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The "improved dynamics" of loop quantum cosmology is extended to include 
anisotropics of the Bianchi I model. As in the isotropic massless scalar 

field serves as a relational time parameter. However, the extension is non-trivial 
because one has to face several conceptual subtleties as well as technical difficulties. 
These include: a better understanding of the relation between loop quantum gravity 
(LQG) and loop quantum cosmology (LQC); handling novel features associated with 
the non-local field strength operator in presence of anisotropics; and finding dynam- 
ical variables that make the action of the Hamiltonian constraint manageable. Our 
analysis provides a conceptually complete description that overcomes limitations of 
earlier works. We again find that the big bang singularity is resolved by quantum 
geometry effects but, because of the presence of Weyl curvature, Planck scale physics 
is now much richer than in the isotropic case. Since the Bianchi I models play a key 
role in the Belinskii, Khalatnikov, Lifshitz (BKL) conjecture on the nature of generic 
space-like singularities in general relativity, the quantum dynamics of Bianchi I cos- 
mologies is likely to provide considerable intuition about the fate of generic space-like 
singularities in quantum gravity. Finally, we show that the quantum dynamics of 
Bianchi I cosmologies projects down exactly to that of the Friedmann model. This 
opens a new avenue to relate more complicated models to simpler ones, thereby 
providing a new tool to relate the quantum dynamics of LQG to that of LQC. 

PACS numbers: 98.80Qc,04.60.Pp, 04.60.-m 



I. INTRODUCTION 

Loop quantum gravity (LQG) [1-3] is a non-perturbative, background independent ap- 
proach to the unification of general relativity and quantum physics. One of its key features is 
that space-time geometry is treated quantum mechanically from the beginning. Loop quan- 
tum cosmology (LQC) [4, 5] is constructed by applying methods of LQG to mini-superspaces 
obtained by a symmetry reduction of general relativity. In the homogeneous, isotropic cos- 
mological models with a massless scalar field, quantum geometry effects of LQG have been 
shown to create a new repulsive force in the Planck regime. The force is so strong that the 
big bang is replaced by a specific type of quantum bounce. The force rises very quickly once 
the scalar curvature reaches ~ — 0.157r/^|j (or matter density p reaches ~ 0.01 ppi) to cause 
the bounce but also dies very quickly after the bounce once the scalar curvature and the 
density fall below these values. Therefore outside the Planck regime the quantum space-time 
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of LQC is very well approximated by the space-time continuum of general relativity. This 
scenario is borne out in the k=0, A=0 models, [6-13], A 7^0 models [14, 15], the k=l closed 
model [16, 17], k=— 1 open model [18] and the k=0 model with an inflationary potential 
with phenomenologically viable parameters [19]. Going beyond the big-bang and big-crunch 
singularities, LQC has also been used to argue that its quantum geometry effects resolve 
all strong curvature singularities in homogeneous, isotropic situations in which matter is 
a perfect fluid with an equation of state of the standard type, p = p{p) [20]. (For recent 
reviews, see, e.g., [21, 22].) Finally, recent investigations [23, 24] of Gowdy models, which 
have an infinite number of degrees of freedom, also indicate that the big-bang is replaced by 
a quantum bounce. 

Detailed and viable quantum theories were constructed in the homogeneous, isotropic 
case using the so-called "/i" scheme. A key open question has been whether or not the 
qualitative features of their Planck scale physics will persist in more realistic situations in 
which these strong symmetry assumptions do not hold exactly. A first step in this direction is 
to retain homogeneity and extend the "improved dynamics" of [10] to anisotropic situations. 
In the isotropic case, there is only one non-trivial curvature invariant, the (space-time) scalar 
curvature (or, equivalently, matter density). In anisotropic situations Weyl curvature is non- 
zero and it too diverges at the big bang. Therefore, now one can enter the Planck regime in 
several inequivalent ways which suggests that the Planck scale physics would now be much 
richer. 

In this paper we will continue the LQC explorations of this issue by analyzing in detail 
the simplest of anisotropic models, the Bianchi I cosmologies. (Previous work on this model 
is discussed below.) As in the isotropic case we will use a massless scalar field as the matter 
source, and it will continue to provide the "relational" or "internal" time a la Leibniz with 
respect to which other physical quantities of interest — e.g., curvatures, shears, expansion 
and matter density — "evolve" . Again, as in the isotropic case, the framework can be further 
extended to accommodate additional matter fields in a rather straightforward fashion. 

Although the Bianchi I models are the simplest among anisotropic cosmologies, results 
obtained in the context of the Belinskii, Khalatnikov, Lifshitz (BKL) conjecture [25, 26] 
suggest that they are perhaps the most interesting ones for the issue of singularity resolu- 
tion. The BKL conjecture states that, as one approaches space-hke singularities in general 
relativity, terms with time derivatives would dominate over those with spatial derivatives, 
implying that the asymptotic dynamics would be well described by an ordinary differential 
equation. By now considerable evidence has accumulated in favor of this conjecture [27-31]. 
For the case when the matter source is a massless scalar field in full general relativity with- 
out any symmetry assumption, these results suggest that, as the system enters the Planck 
regime, dynamics along any fixed spatial point would be well described by a Bianchi I met- 
ric. Therefore understanding the fate of Bianchi I models in LQC could provide substantial 
intuition for what happens to generic space- like singularities in LQG [32, 33]. 

Indeed, in cosmological contexts where one has approximate homogeneity, a natural strat- 
egy in full LQG is to divide the spatial 3-manifold into small, elementary cells and assume 
that there is homogeneity in each cell, with fields changing slowly as one moves from one 
cell to the next. (For an exploration along these lines in the older "/^o scheme," see [34].) 
Now, if one were to assume that geometry in each elementary cell is also isotropic, then the 
Weyl tensor in each cell — and therefore everywhere — would be forced to be zero. A natural 
strategy to accommodate realistic, non-vanishing Weyl curvature would be to use Bianchi 
I geometry in each cell and let the parameters ki vary slowly from one cell to another. In 
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this manner, LQC of the Bianchi I model can pave way to the analysis of the fate of generic 
space-like singularities of general relativity in full LQG. 

Because of these potential applications, Bianchi I models have already drawn considerable 
attention in LQC (see in particular [35-40]). During these investigations, groundwork was 
laid down which we will use extensively. However, in the spatially non-compact context 
(i.e., when the spatial topology is M.^ rather than T^), the construction of the quantum 
Hamiltonian constraint turned out to be problematic. The Hamiltonian constraint used 
in the early work has the same difficulties as those encountered in the "//o-scheme" in the 
isotropic case (see, e.g., [12], or Appendix B of [21]). More recent papers have tried to 
overcome these limitations by mimicking the scheme used successfully in the isotropic 
case. However, to make concrete progress, at a key point in the analysis a simplifying 
assumption was made without a systematic justification.^ Unfortunately, it leads to quantum 
dynamics which depends, even to leading order, on the choice of an auxiliary structure (i.e., 
the fiducial cell) used in the construction of the Hamiltonian framework [40]. This is a 
major conceptual drawback. Also, the final results inherit certain features that arc not 
physically viable (e.g. the dependence of the quantum bounce on "directional densities" in 
[36, 37]). We will provide a systematic treatment of quantum dynamics that is free from 
these drawbacks. 

To achieve this goal one has to overcome rather non-trivial obstacles which had stalled 
progress for the past two years. This requires significant new inputs. The first is conceptual: 
we will sharpen the correspondence between LQG and LQC that underlies the definition 
of the curvature operator F*^ in terms of holonomies. The holonomies we are led to use in 
this construction will have a non-trivial dependence on triads, stemming from the choice 
of loops on which they are evaluated (see footnote 1). As a result, at first it seems very 
difficult to define the action of the resulting quantum holonomy operators. Indeed this 
was the primary technical obstacle that forced earlier investigations to take certain short 
cuts — the assumption mentioned above — while defining F*^. The second new input is the 
definition of these holonomy operators without having to take a recourse to such short cuts. 
But then the resulting Hamiltonian constraint appears unwieldy at first. The third major 
input is a rearrangement of configuration variables that makes the constraint tractable both 
analytically, as in this paper, and for the numerical work in progress [41]. 

Finally, we will find that the resulting Hamiltonian constraint has a striking feature which 
could provide a powerful new tool in relating the quantum dynamics of more complicated 
models to that of simpler models. It turns out that, in LQC, there is a well-defined projection 
from the Bianchi I physical states to the Priedmann physical states which maps the Bianchi 
I quantum dynamics exactly to the isotropic quantum dynamics. Previous investigations of 
the relation between quantum dynamics of a more complicated model to that of a simpler 
model generally began with an embedding of the Hilbert space T^Res of the more restricted 
model in the Hilbert space Hcen of the more general model (see, e.g., [43, 44]). In generic 
situations, the image of T^Res under this embedding was not left invariant by the more general 
dynamics on Tioen- This led to a concern that the physics resulting from first reducing and 



In the isotropic case, "improved" dynamics [10] required that fi be proportional to l/-\/[p|. In the 
anisotropic case, one has three pi and quantum dynamics requires the introduction of three fii. In the 
Bianchi I case now under consideration, it was simply assumed [36, 37, 40] that jli be proportional to 
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then quantizing may be completely different from that obtained by quantizing the larger 
system and regarding the smaller system as its sub-system. The new idea of projecting from 
'Hcen to T^Res corrcspouds to "integrating out the degrees of freedom that are inaccessible 
to the restricted model" while the embedding Hrss in to Ti-oen corresponds to "freezing by 
hand" these extra degrees of freedom. Classically, both are equally good procedures and 
in fact the embedding is generally easier to construct. However, in quantum mechanics it 
is more appropriate to integrate out the "extra" degrees of freedom. In the present case, 
one "integrates out" anisotropies to go from the LQC of the Bianchi I models to that of 
the Priedmann model. This idea was already proposed and used in [42] in a perturbative 
treatment of anisotropies in locally rotationally symmetric, diagonal, Bianchi I model. We 
extend that work in that we consider the full quantum dynamics of diagonal Bianchi I model 
without additional symmetries and, furthermore, use the analog of the "/i scheme" in which 
the quantum constraint is considerably more involved than in the "/^o-type" scheme used in 
[42]. The fact that the LQC dynamics of the Friedmann model is recovered exactly provides 
some concrete support for the hope that LQC may capture the essential features of full 
LQG, as far as the quantum dynamics of the homogeneous, isotropic degree of freedom is 
concerned. 

The material is organized as follows. We will begin in section II with an outline of the 

classical dynamics of Bianchi I models. This overview will not be comprehensive as our goal 
is only to set the stage for the quantum theory which is developed in section 111. In section 
IV we discuss three key properties of quantum dynamics: the projection map mentioned 
above, agreement of the LQC dynamics with that of the Wheeler DeWitt theory away 
from the Planck regime and effective equations. (The isotropic analogs of these equations 
approximate the full LQC dynamics of Friedmann models extremely well.) In section IV 
we summarize the main results and discuss some of their ramifications. The Appendix A 
discusses parity type discrete symmetries which play an important role in the analysis of 
quantum dynamics. 

II. HAMILTONIAN FRAMEWORK 

In this section we will summarize those aspects of the classical theory that will be needed 
for quantization. For a more complete description of the classical dynamics see, e.g., [35- 
37, 45]. 

Our space-time manifold M will be topologically R^. As is standard in the literature on 
Bianchi models, we will restrict ourselves to diagonal Bianchi I metrics. Then one can fix 
Cartesian coordinates r, Xi on M and express the space-time metric as: 

ds'^ = -N^dr'^ + a\ dx\ + aj dxj + dxj , (2.1) 

where is the lapse and are the directional scale factors. Thus, the dynamical degrees 
of freedom are encoded in three functions aj(T) of time. Bianchi I symmetries permit us 
to rescale the three spatial coordinates Xi by independent constants. Under the 
directional scale factors transform as^ Oj a^^Ui. Thus, the numerical value of a directional 
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scale factor, say ai, is not an observable; only ratios such as ai(r)/ai(T') are. The matter 
source will be a massless scalar field which will serve as the relational or internal time. 
Therefore, it is convenient to work with a harmonic time function, i.e. to ask that r satisfy 
□r = 0. Prom now on we will work with this choice. 

Since the spatial manifold is non-compact and all fields are spatially homogeneous, to 
construct a Lagrangian or a Hamiltonian framework one has to introduce an elementary cell 
V and restrict all integrations to it [7] . We will choose V so that its edges lie along the fixed 
coordinate axis Xi. As in the isotropic case, it is also convenient to fix a fiducial fiat metric 
"Qab with line element 

dsl = dxj + dxl + dxl . (2.2) 

We will denote by °q the determinant of this metric, by Li the lengths of the three edges of V 
as measured by "qah^ and by Vo = L1L2L3 the volume of the elementary cell V also measured 
using "Qab- Finally, we introduce fiducial co-triads lu^ = DaX"^ and the triads "e" dual to 
them. Clearly they are adapted to the edges of V and are compatible with °qab (i.e., satisfy 
°Qab — As noted above, Bianchi I symmetries allow each of the three coordinates 

to be rescaled by an independent constant ai. Under these rcscalings, 
co-triads transform as = ai°u!l, and triads °e" are rescaled by inverse powers of a,. The 
fiducial metric is transformed to "g^^ defined by ds'^ := a^dxf + a^dx"^ + aldx\. We must 
ensure that our physical results do not change under these rescalings. Finally, the physical 
co-triads are given by cj* = a*°cj* and the physical 3- metric gab is given by Qab — ^ll^l ^ij- 

With these fiducial structures at hand, we can now introduce the phase space. Recall 
first that in LQG the canonical pair consists of an SU(2) connection and a triad Ef of 
density weight one. Using the Bianchi I symmetry, from each gauge equivalence class of 
these pairs we can select one and only one, given by: 

A^, =: diVy' and ^ ^ge" =: p,L,V-'y^%t, (2-3) 

where Ci,pi are constants and q = (P1P2P3) °g KT^ is the determinant of the physical spatial 
metric Qab- Thus the connections A\ are now labelled by three parameters d and the triads 
Ef by three parameters pi. If pi are positive, the physical triad ef and the fiducial triad 
°ef have the same orientation. A change in sign of, say, pi corresponds to a change in the 
orientation of the physical triad brought about by the flip e" —e^. These flips are gauge 
transformations because they do not change the physical metric Qab- The momenta Pi are 
directly related to the directional scale factors: 

pi = sgn(ai)|a2a3| L2L3, p2 = sgn(a2)|aia3| L1L3, = sgn(a3)|aia2| L1L2 . (2.4) 

where we take the directional scale factor to be positive if the triad vector is parallel 
to °e" and negative if it is anti-parallel. As we will see below, in any solution to the field 
equations, the connection components q are directly related to the time derivatives of Oj. 

The factors of Li in (2.3) ensure that this parametrization is unchanged if the fiducial 
co-triad, triad and metric are rescaled via Xi aiXi. However, the parametrization does 
depend on the choice of the cell V. Thus the situation is the same as in the isotropic case 
[7]. (The physical fields and Ef are of course insensitive to changes in the fiducial metric 
or the cell.) To evaluate the symplectic structure of the symmetry reduced theory, as in the 
isotropic case [7], we begin with the expression of the symplectic structure in the full theory 
and simply restrict the integration to the cell V. The resulting (non-vanishing) Poisson 
brackets are given by: 

{c^p,■} = 877^7 5}. (2.5) 
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To summarize, the phase space in the Bianchi I model is six dimensional, coordinatized 
by pairs c*,pj, subject to the Poisson bracket relations (2.5). This description is tied to the 
choice of the fiducial cell V but is insensitive to the choice of fiducial triads, co-triads and 
metrics. 

Next, let us consider constraints. The full theory has a set of three constraints: the 
Gauss, the diffeomorphism and the Hamiltonian constraints. It is straightforward to check 
that, because we have restricted ourselves to diagonal metrics and fixed the internal gauge, 
the Gauss and the diffeomorphism constraints are identically satisfied. We are thus left with 
just the Hamiltonian constraint. Its expression is obtained by restricting the integration in 
the full theory to the fiducial cell V: 

Ch — Cgrav + Cmatt = N (Hgrav + 'Hmatt) d^^; (2.6) 

Jv 

where N is the lapse function and the gravitational and the matter parts of the constraint 
densities are given by 



nj * J 

'l-grav — 



167rG'. 



(e'hFj - 2(1 + j^)KiK0 and Ti^att = Vq Pmatt- (2.7) 



Here 7 is the Barbero-Immirzi parameter, F^b'' is the curvature of the connection A^, given 

by 

F„b'= = 29[„A,]^ + e,/A^A^,, (2.8) 

i^* is related to the extrinsic curvature Kab via fC* = Kab^^^ and Pmatt is the energy density 
of the matter fields. In general, is related to K^^ and the spin connection defined by 
the triad ef via = T^ + '^K^. However, because Bianchi I models are spatially flat, = 
in the gauge chosen in (2.3), whence A\ — jK^. This property and the fact that spatial 
derivatives of vanish by the Bianchi I symmetry leads us to the relation 

= r'e'^kFab' . (2.9) 

Therefore, the gravitational part of the Hamiltonian constraint can be simplified: 



'ij ' ■/ 'ij rp k 



(P1P2C1C2 + P1P3C1C3 + P2P3C2C3) , (2-10) 



Finally, recall that our matter field is a massless scalar field T. The matter energy density 
of the scalar field T is given by pmatt = P(t)/2V^^, where V = ^/ipilhPsl is the physical volume 
of the elementary cell. Our choice of harmonic time r implies that the lapse function is given 
by = ^/\plP2P3\■ With these choices the constraint (2.6) simplifies further: 



Ch= - .r^:" e.,'F^ + d^x (2.11) 




1 

(P1P2C1C2 + P1P3C1C3 + P2P3C2C3) + (2-12) 
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Physical states of the classical theory lie on the constraint surface Ch = 0. The time 
evolution of each pi and q is obtained by taking their Poisson bracket with Ch- 

^ = {p„Ch} = -SnGj^ = ^ (P2C2 +P3C3) ; (2.13) 

QT OCi 7 

^ = {ci,C^,} = SttGt^ = — (P2C2 +P3C3) . (2.14) 
dr api 7 

The four other time derivatives can be obtained via permutations. Although the phase 
space coordinates d,Pi themselves depend on the choice of the fiducial cell V, the dynamical 
equations for and E'^ — and hence also for the physical metric Qah and the extrinsic 
curvature Kab — that follow from (2.13) and (2.14) are independent of this choice. 
Combining Eqs. (2.4), (2.13) and (2.14), one finds 

a = 7L, (010203)-' ^ (2.15) 

dr 

It is instructive to relate the Cj to the directional Hubble parameters Hi = dlnoj/dt where 
t is the proper time, corresponding to the lapse function N(^t) = 1- Since t is related to the 
harmonic time r via Ndr = iV(t)dt 

(2.16) 



Therefore, we have 

Q = 7Li--^ = -fLiUiHi (2.17) 

where LjOj is the length of the ith edge of V as measured by the physical metric Qab- 

Next, it is convenient to introduce a mean scale factor a :— (010203) ^''^ which encodes 
the physical volume element but ignores anisotropies. Then, the mean Hubble parameter is 
given by 

H :=^^ = 1 (H1 + H2 + H3), where as before Hr.= ^^. (2.18) 
at 3 at 

Squaring Eq. (2.18) and using the implication 

H,H2 + H^Hs + HsH^ = Stt Gp^^tt (2.19) 

of the Hamiltonian constraint, we obtain the generalized Priedmann equation for Bianchi I 
space-times, 

H'^—-Pm.u + —, (2.20) 

where 

= ^ [{Hr - H,f + {H, - + {H^ - H,)'] (2.21) 

is the shear term. The right hand side of (2.20) brings out the fact that the anisotropic 
shears {Hi — Hj) contribute to the energy density; they quantify the energy density in the 
gravitational waves. Using the fact that our matter field has zero anisotropic stress one 
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can show that T? is a constant of the motion [37]. If the space-time itself is isotropic, 
then = and Eq. (2.20) reduces to the usual Friedmann equation for the standard 
isotropic cosmology. These considerations will be useful in interpreting quantum dynamics 
and exploring the relation between the Bianchi I and Friedmann quantum Hamiltonian 
constraints. 

Next, let us consider the scalar field T. Because there is no potential for it, its canonically 
conjugate momentum p{T) is a constant of motion (which, for definiteness, will be assumed 
to be positive). Therefore, in any solution to the field equations T grows linearly in the 
harmonic time r. Thus, although T does not have the physical dimensions of time, it is 
a good evolution parameter in the classical theory. The form of the quantum Hamiltonian 
constraint is such that T will also serve as a viable internal time parameter in the quantum 
theory. 

We will conclude with a discussion of discrete 'refiection symmetries' that will play an 
important role in the quantum theory. (For further details see the Appendix.) In the 
isotropic case, there is a single refiection symmetry, 11 (p) = —p which physically corresponds 
to the orientation reversal e" — > — of triads. These are large gauge transformations, 
under which the metric qah remains unchanged. The Hamiltonian constraint is invariant 
under this reflection whence one can, if one so wishes, restrict one's attention just to the 
sector p > of the phase space. In the Bianchi I case, we have three reflections Hj, each 
corresponding to the flip of one of the triad vectors, leaving the other two untouched (e.g., 
ni(pi,P2,P3) = {—Pi-iPi-iP^))- As shown in [46], the Hamiltonian flow is left invariant under 
the action of each H,. Therefore, it suffices to restrict one's attention to the positive octant 
in which all three Pi are non-negative: dynamics in any of the other seven octants can be 
easily recovered from that in the positive octant by the action of the discrete symmetries 
H,. 

• Remark: In the LQC literature on Bianchi I models, a physical distinction has occa- 
sionally been made between the fiducial cells V which are "cubical" with respect to 
the fiducial metric °qab and those that are "rectangular." (In the former case all Lj are 
equal.) However, given any cell V one can always find a flat metric in our collection 
(2.1) with respect to which that V cubical. Using it as °qah one would be led to call 
it cubical. Therefore the distinction is unphysical and the hope that the restriction to 
cubical cells may resolve some of the physical problems faced in [36, 37] was misplaced. 

III. QUANTUM THEORY 

This section is divided into four parts. In the first, we briefiy recall quantum kinemat- 
ics, emphasizing issues that have not been discussed in the literature. In the second, we 
spell out a simple but well-motivated correspondence between the LQG and LQC quantum 
states that plays an important role in the definition of the curvature operator Fa^^ in terms 
of holonomies. However, the paths along which holonomies are evaluated depend in a rather 
complicated way on the triad (or momentum) operators, whence at first it seems very dif- 
ficult to define these holonomy operators. In the third subsection we show that geometric 
considerations provide a natural avenue to overcome these apparent obstacles. The resulting 
Hamiltonian constraint is, however, rather unwieldy to work with. In the last sub-section 
we make a convenient redefinition of configuration variables to simplify its action. The sim- 
plification, in turn, will provide the precise sense in which the singularity is resolved in the 
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quantum theory. 



A. LQC Kinematics 



We will summarize quantum kinematics only briefly; for details, see e.g. [36, 37]. Let 
us begin by specifying the elementary functions on the classical phase space which are to 
have unambiguous analogs in the quantum theory. In LQC this choice is directly motivated 
by the structure of full LQG [1-3]. As one might expect from the isotropic case [7, 9], the 
elementary variables are the three momenta pi and holonomics /i^ along edges parallel to 
the three axis Xi, where iLi is the length of the edge with respect to the fiducial metric °qab-^ 
These functions are (over) complete in the sense that they suffice to separate points of the 
phase space. Taking the xi axis for concreteness, the holonomy hf^ has the form 

h\'{ci,C2,C3) = COS — I + 2 sm — Ti (3.1) 

where I is the unit 2x2 matrix and Tj constitute a basis of the Lie algebra of SU(2), satisfying 
rV-' = \e^^ — \5'^^1. Thus, the holonomies are completely determined by almost periodic 
functions exp(i£cj) of the connection; they are called "almost" periodic because i is any real 
number rather than an integer. In quantum theory, then, elementary operators hf' and pi 
are well-defined and our task is to express other operators of physical interest in terms of 
these elementary ones. 

Recall that in the isotropic case it is simplest to specify the gravitational sector of the 
kinematic Hilbert space in the triad of p representation: it consists of wave functions \£'(p) 
which arc symmetric under p —p and have a finite norm: = |\E'(p)p < oo. In 

the Bianchi I case it is again simplest to describe H.^^ in the momentum representation. 
Consider first a countable linear combination, 

l^)= *(Pl,P2,P3)|Pl,P2,P3) with |^'(pi,p2,P3)r < OO, (3.2) 

Pl,P2,P3 Pl,P2,P3 

of orthonormal basis states \pi,P2,'P3)^ where 

(Pl,P2,P3 |pi,P2,P3) = ^Pip'i ^P2Pi ^P3P^ ■ (3-3) 

Next, recall that on the classical phase space the three refiections 11 j represent large gauge 
transformations under which physics does not change. They have a natural induced ac- 
tion Ilj on the space of wave functions ^'(pi,P2,P3)- (Thus, for example, IIi^(pi,p2,P3) = 
^(— Pi,P2,P3)-) Physical observables commute with flj. Therefore, as in gauge theories, 
each eigenspace of ftj provides a physical sector of the theory. Since ft? = I, eigenvalues of 
Ilj are ±1. For definiteness, as in the isotropic case, we will assume that the wave functions 
^{PiiP2)P^) arc symmetric under Ilj. Thus, the gravitational part "H^^ of the kinematical 
Hilbert space is spanned by wave functions \E'(pi,p2,P3) satisfying 

*(Pl,P2,P3) = *(|pi|, \P2\, \P2\) (3-4) 



^ More precisely, the dimensionlcss number I is the length of the edge along which the holonomy is evaluated, 
measured in the units of the length of the edge of V parallel to it. Since £ is a ratio of lengths, its value 
does not depend on the fiducial or any other metric. 
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which have finite norm (3.2). 

The basis states \pi,P2,P3) are eigenstates of quantum geometry: In the state \pi,P2,P3) 
the face Si of the fiducial cell V orthogonal to the axis Xi has area \pi \ . Note that although 
Pi e ]R, the orthonormality holds via Kronecker deltas rather than the usual Dirac distribu- 
tions; this is why the LQC quantum kinematics is inequivalent to that of the Schrodinger 
theory used in Wheeler DeWitt cosmology. Finally the action of the elementary operators 
is given by: 

Pi \Pl:P2:P3) = Pi \Pl:P2:P3) and CXp i^Ci P2, Ps) = \Pl " Sn-fGM, P2,P3) (3.5) 

and similarly for p2 , exp iic2 , Ps and exp iic-^ . 

The full kinematical Hilbert space TYkin will be the tensor product, TCkin = ^lin^ ® ^Sn" 
where, as in the isotropic case, we will set '^^Sn" = L^(R, dT) for the Hilbert space of 
the homogeneous scalar field T. On H^^^, the operator T will act by multiphcation and 
P{T) '■= —ihd/dT will act by differentiation. Note that we can also use a "polymer Hilbert 
space" for spanned by almost periodic functions of T. The quantum Hamiltonian 

constraint (3.22) will remain unchanged and our construction of the physical Hilbert space 
will go through as it is [47] . 

B. The curvature operator Fab'' 

To discuss quantum dynamics, we have to construct the quantum analog of the Hamil- 
tonian constraint. Since there is no operator corresponding to the connection coefficients q 
on TYkin^' cannot use (2.12) directly. Rather, as in the isotropic case [10], we will return 
to the expression (2.11) involving curvature Fab'^. Our task then is to find the operator on 
^kin^ corresponding to Fab''- As is usual in LQG, the idea is to first express the curvature 
in terms of our elementary variables — holonomies and triads — and then replace them by 
their direct quantum analogs. Recall first that, in the classical theory, the a-b component 
of Fab'' can be written in terms of holonomies around a plaquette (i.e., a rectangular closed 
loop whose edges are parallel to two of the axes Xi): 

Fab" ^2 lim Tr 

where Ara is the area of the plaquette □ and the holonomy hoij around the plaquette Dij 
is given by 

=/,f )"\f^)-\f (3.7) 

where fij Lj is the length of the jth edge of the plaquette, as measured by the fiducial metric 
°qab- (There is no summation over i,j.) Because the Ara is shrunk to zero, the limit is 
not sensitive to the precise choice of the closed plaquette □. Now, in LQG the connection 
operator does not exist, whence if we regard the right side of (3.6) as an operator, the 
limit fails to converge in H.^^^ . The non-existence of the connection operator is a direct 
consequence of the underlying diffeomorphism invariance [48] and is intertwined with the fact 
that the eigenvalues of geometric operators — such as the area operator Arc\ associated with 
the plaquette under consideration — are purely discrete. Therefore, in LQC the viewpoint is 
that the non-existence of the limit Arn — > in quantum theory is not accidental: quantum 
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(3.6) 
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(a) (b) 

FIG. 1: Depiction of the LQG quantum geometry state corresponding to the LQC state \pi,P2,P3)- 
The LQG spin-network has edges parallel to the three axes selected by the diagonal Bianchi I 
symmetries, each carrying a spin label j = 1/2. (a) Edges of the spin network traversing through 
the fiducial cell V. (b) Edges of the spin network traversing the 1-2 face of V and an elementary 
plaquette associated with a single flux line. This plaquette encloses the smallest quantum, A^pj, 
of area. The curvature operator ^12*^ is defined by the holonomy around such a plaquette. 



geometry is simply telling us that we should shrink the plaquette not till the area it encloses 
goes to zero, but rather only to the minimum non-zero eigenvalue A fpj of the area operator 
(where A is a dimensionless number). The resulting quantum operator Fab^ then inherits 
Planck scale non-localities. 

To implement this strategy in full LQG one must resolve a difficult issue. If the plaquette 
is to be shrunk only to a finite size, the operator on the right side of (3.6) would depend 
on what that limiting plaquette is. So, which of the many plaquettes enclosing an area 
A£|[ should one use? Without a well-controlled gauge fixing procedure, it would be very 
difficult to single out such plaquettes, one for each 2-dimensional plane in the tangent space 
at each spatial point. However, in the diagonal Bianchi I case now under consideration, a 
natural gauge fixing is available and indeed we have already carried it out. Thus, in the 
i-j plane, it is natural to choose a plaquette Dij so that its edges are parallel to the Xi-Xj 
axis. Furthermore, the underlying homogeneity implies that it suffices to introduce the three 
plaquettes at any one point in our spatial 3-manifold. 

These considerations severely limit the choice of plaquettes but they do not determine 
the lengths of the two edges in each of these plaquettes. To completely determine the 
plaquettes, as in the isotropic case, we will use a simple but well-motivated correspondence 
between kinematic states in LQG and those in LQC. However, because of anisotropies, 
new complications arise which require that the correspondence be made much more precise. 
Fix a state \pi,P2,P3) in '^lin^ LQC. In this state, the three faces of the fiducial cell 
V orthogonal to the Xj-axis have areas \pi\ in the LQC quantum geometry. This is the 
complete physical information in the ket \pi,P2,P3)- How would this quantum geometry be 
represented in full LQG? First, the macroscopic geometry must be spatially homogeneous 
and we have singled out three axes with respect to which our metrics are diagonal. Therefore, 
semi-heuristic considerations suggest that the corresponding LQG quantum geometry state 
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should be represented by a spin network consisting of edges parallel to the three axes (see 
Fig. 1(a)). Microscopically this state is not exactly homogeneous. But the coarse grained 
geometry should be homogeneous. To achieve the best possible coarse grained homogeneity, 
the edges should be packed as tightly as is possible in the desired quantum geometry. That 
is, each edge should carry the smallest non-zero label possible, namely j = 1/2. 

For definiteness, let us consider the 1-2 face Su of the fiducial cell V which is orthogonal 
to the xs axis (see Fig. 1(b)). Quantum geometry of LQG tells us that at each intersection 
of any one of its edges with 5'i2, the spin network contributes a quantum of area A ip^ on this 
surface, where A = 47r7-\/3 [49]. For this LQG state to reproduce the LQC state |pi,P2,P3) 
under consideration 5'i2 must be pierced by A^3 edges of the LQG spin network, where A^3 is 
given by 

NsA£l,= \p,\. 

Thus, we can divide 5*12 into identical rectangles each of which is pierced by exactly one 
edge of the LQG state, as in Fig. 1(b). Any one of these elementary rectangles encloses 
an area A^lj and provides us the required plaquette ^12. Let the dimensionless lengths of 
the edges of these plaquettes be fii and jl2- Then their lengths with respect to the fiducial 
metric °qab are fiiLi and /i2-f'2- Since the area of 5'i2 with respect to °qab is L1L2, we have 

A^3 jliLijl2L2 = L1L2 . 

Equating the expressions of A^a from the last two equations, we obtain 

^l^^l2 = ^ . (3.8) 

This relation by itself docs not fix /2i and fi2- However, repeating this procedure for the 2-3 
face and the 3-1 face, we obtain, in addition, two cyclic permutations of this last equation 
and the three simultaneous equations do suffice to determine fif 



- _ bi|A4i - _ b2|A4, _ balA^i 

V IP2P3I V \pm\ V IP1P2I 

To summarize, by exploiting the Bianchi I symmetries and using a simple but well- 
motivated correspondence between LQG and LQC states we have determined the required 
elementary plaquettes enclosing an area A £pj on each of the three faces of the cell V. On 
the face Sij, the plaquette is a rectangle whose sides are parallel to the Xi and X j axes 
and whose dimensionless lengths are pi and flj respectively, given by (3.9). Note that (as 
in the isotropic case [10]) the /2j and hence the plaquettes are not fixed once and for all; 
they depend on the LQC state \pi,P2,P3) of quantum geometry in a specific fashion. The 
functional form of this dependence is crucial to ensure that the resulting quantum dynamics 
is free from the difficulties encountered in earlier works. 

Components of the curvature operator F^b'' can now be expressed in terms of holonomies 
around these plaquettes: 

Fab' = 2 Tr (^^^^ r'^ X (3-10) 

with 

/in . = hf^^~' hf'^~' hf^ hf'^ , (3.11) 
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where Hj are given by (3.9). (There is no summation over Using the expression (3.1) 

of holonomies, it is straightforward to evaluate the right hand side. One finds: 

F,'=,./(!!^».„y(!i^«.,y, (3.12) 

where the usual summation convention for repeated covariant and contravariant indices 
applies and 



where there is now no sum over i. This is the curvature operator we were seeking. 

We will conclude with a discussion of the important features of this procedure and of the 
resulting quantum dynamics. 

1. In the isotropic case all pi are equal {pi — p) whence our expressions for /ij reduce 
to a single formula, Jl = ^J'K^J]p\. This is precisely the result that was obtained in the 
"improved dynamics" scheme for the k = isotropic models. Thus, we have obtained a 
generalization of that result to Bianchi I models. 

2. In both cases, the key observation is that the plaquette should be shrunk till its area 
with respect to the physical — rather than the fiducial — geometry is A^pj. However, there 
are also some differences. First, in the above analysis we set up and used a correspondence 
between quantum geometries of LQG and LQC in the context of Bianchi I models. In 
contrast to the previous treatment in the isotropic models [10], we did not have to bring in 
classical geometry in the intermediate steps. In this sense, even for the isotropic case, the 
current analysis is an improvement over what is available in the literature. 

3. A second difference between our present analysis and that of [10] is the following. Here, 
the semi-heuristic representation of LQC states \pi,P2,P3) in terms of spin networks of LQG 
suggested that we should consider spin networks which pierce the faces of the fiducial cell 
V as in Fig. 1(a). (As one would expect, these states are gauge invariant.) The minimum 
non-zero eigenvalue of the area operator on such states is Aip^ with A = A^/Sn^f. This is 
twice the absolute minimum of non-zero eigenvalue on all gauge invariant states. However, 
that lower value is achieved on spin networks (whose edges are again labelled hy j — 1/2 but) 
which do not pierce the surface but rather intersect it from only one side. (In order for the 
state to be gauge invariant, the edge then has to continue along a direction tangential to the 
surface. For details, see [49].) Obvious considerations suggest that such states cannot feature 
in homogeneous models. Since the discussion in the isotropic case invoked a correspondence 
between LQG and LQC at a rougher level, this point was not noticed and the value of A 
used in [10] was 2\^-k^. We emphasize, however, that although the current discussion is 
more refined, it is not a self-contained derivation. A more complete analysis may well change 
this numerical factor again. 

4. On the other hand, we believe that the functional dependence of fii on pi is robust: 
As in the isotropic case this dependence appears to be essential to make quantum dynamics 
viable. Otherwise quantum dynamics can either depend on the choice of the fiducial cell V 
even to leading order, or is physically incorrect because it allows quantum effects to dominate 
in otherwise "tame" situations, or both. The previous detailed, quantum treatments of the 
Bianchi I model in LQC did not have this functional dependence because they lacked the 
correspondence between LQG and LQC we used. Rather, they proceeded by analogy. As 



14 



we noted above, in the isotropic case there is a single n and a single p and the two are 
related hy p, = A/Afpj/|p|. The most straightforward generalization of this relation to 
Bianchi I models is fii = \J A£pj/ \pi \ . This expression was simply postulated and then used 
to construct quantum dynamics [36, 37]. The resulting analysis has provided a number 
of useful technical insights. However, this quantum dynamics suffers from the problems 
mentioned above [40] . The possibility that the correct generalization of the isotropic results 
to Bianchi I models may be given by (3.9) was noted in [38, 50] and in the Appendix C of 
[37]. However, for reasons explained in the next sub-section, construction of the quantum 
Hamiltonian operator based on (3.9) was thought not to be feasible. Therefore, this avenue 
was used only to gain qualitative insights and was not pursued in the full quantum theory. 



C. The quantum Hamiltonian constraint 

With the curvature operator Fab' at hand, it is straightforward to construct the quan- 
tum analog of the Hamiltonian constraint (2.6) because the triad operators can be readily 
constructed from the three pi. Ignoring for a moment the factor-ordering issues, the gravi- 
tational part of this operator is given by 



Cgrav = - o-./^„?2 A p2 \PiP'i\P^\ A«lCl Slu //2C2 + Pl\p2\Pz Slu //iCi Siu /igCa 



+ IP2P3 sin //2C2 sin ^303] 



(3.14) 



where for simplicity of notation here and in what follows we have dropped hats on p^ and 
sin/iiQ. To write the action of this operator on V-^^ , it suffices to specify the action of the 
operators exp(i/ijCi) on the kinematical states ^(pi,p2,P3)- The expression (3.9) of fxi and 
the Poisson brackets (2.5) imply: 



exp(±i//iCi) = exp ^ ^ 



Pi 



P2P5 



d 

dpi 



(3.15) 



and its cyclic permutations. At first sight this expression seems too complicated to yield a 
manageable Hamiltonian constraint. 



• Remark: In the isotropic case, the corresponding expression is simply 



exp(±i//c) = exp ^ =p 8777 

Since ^ ~ where v ~ \p\^^'^ is the physical volume of the fiducial cell V, 

this operator can be essentially written as exp(d/dv) and acts just as a displacement 
operator on functions \&(f) of v. In the operator (3.15) by contrast, all three pi 
feature in the exponent. This is why its action was deemed unmanageable. As we 
noted at the end of section IIIB, progress was made [36, 37] by simply postulating an 
alternative, more manageable expression /i, = (\/A ^pi/ \/|Pi|), the obvious analog of 
p, — (s/A £p\) / ^y\p\ in the isotropic case [10]. Then each exp(±i//jCi) can be expressed 
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essentially as a displacement operator expd/d?;j with Vi ~ \pi\^^'^ and the procedure 
used in the isotropic case could be implemented on states '^{vi,V2,V3). Bianchi I 
quantum dynamics then resembled three copies of the isotropic dynamics. However, 
as noted above this solution is not viable [40] . 

Our new observation is that the operator (3.15) can in fact be handled in a manageable 
fashion. Let us first make an algebraic simplification by introducing new dimensionless 
variables Aj : 

(47r|7| v^4i)'/' 

(so that sgn(Ai) = sgn(pj)). Then, we can introduce a new orthonormal basis | Ai, A2, A3) in 
^kin^ by an obvious rescaling. These vectors are again eigenvectors of the operators pi : 

P«|Ai,A2,A3) = sgn(A,)(47r|7|VA4i)iA2|Ai,A2,A3). (3.17) 
We can expand out any ket |^) in H^^^ as j^^) = ^(Ai, A2, A3) |Ai, A2, A3) and re-express the 

— * 

right side of (3.15) as an operator on wave functions ^(A), 

eM±^fi.c,) = exp (^f^ ^) Ef , (3.18) 

where the notation Ef has been introduced as shorthand. (Here, we have used the property 
7 = sgn(pi]92P3)|7| of the Barbero-Immirzi parameter from Appendix A.) To obtain the ex- 

I — * 

plicit action of Ef on wave functions ^(A) we note that, since the operator is an exponential 

— * 

of a vector field, its action is simply to drag the wave function ^(A) a unit affine parameter 
along its integral curves. Furthermore, since the vector field d/dAi is in the Ai direction, 
the coefficient I/A2A3 is constant along each of its integral curves. Therefore it is possible 
to write down the explicit expression of E^ : 

(^f ^) (Ai, A2, A3) = ^(Ai ± A2, A3) . (3.19) 

The non-triviality of this action lies in the fact that while the wave function is dragged 
along the Ai direction, the affine distance involved in this dragging depends on A2, A3. This 
operator is well-defined because our states have support only on a countable number of Aj. 
In particular, the image (^Ef ^f) (A) vanishes identically at points A2 = or A3 = because 
^ does not have support at Ai = 00. Thus the factor A2A3 appearing in the denominator 
does not cause difficulties. 

We can now write out the gravitational part of the Hamiltonian constraint: 

Cgrav — Cgrlv + C^lv + ^graV (3.20) 

with 

^^Iv = -7rMpiA/|AiA2A3| sin p,2C2 sgnA2 | A1A2A3I sgnA3 sin flsCa 

-Fsin/i3C3SgnA3 IA1A2A3I sgnA2sin/i2C2 V^XiA^A^ (3-21) 
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where we have used the simplest symmetric factor ordering that reduces to the one used in 
[11] in the isotropic case. ( CgrL and CgrL are given by the obvious cychc permutations.) In 
Appendix A, we show that, under the action of reflections fij on 7i|[^^, the operators sin /ijCj 
have the same transformation properties that q have under reflections Ilj in the classical 
theory. As a consequence, Cgrav is also reflection symmetric. Therefore, its action is well 
defined on H^^^: Cgrav is a densely defined, symmetric operator on this Hilbert space. In the 
isotropic case, its analog has been shown to be essentially self-adjoint [52]. In what follows 
we will assume that (3.20) is essentially self-adjoint on H^i^ and work with its self-adjoint 
extension. 

Finally, it is straightforward to write down the quantum analog of the full Hamiltonian 
constraint (2.6): 

-h^d^ *(A, T) = e *(A, T) (3.22) 

where = — Cgiav As in the isotropic case, one can obtain the physical Hilbert space Tiphy 
by a group averaging procedure and the result is completely analogous. Elements of Hphy 
consist of 'positive frequency' solutions to (3.22), i.e., solutions to 

-ihdT-^{X,T) ^ ^/\e\■^{X,T) , (3.23) 

which are symmetric under the three refiection maps Ili, i.e. satisfy 

^(Ai, A2, A3, T) = ^(|Ai|, IA2I, IA3I, T) . (3.24) 

The scalar product is given simply by: 

= (*i(A,r„)|*2(A,r„))kin 

= J2 ^i{X,To)M\To) (3.25) 

Ai,A2,A3 

where To is any "instant" of internal time T. 

• Remark: In the isotropic LQC literature [10, 16, 17] one began in the classical theory 
with proper time t (which corresponds to the lapse function N(^t) — 1) and made a 
transition to the relational time provided by the scalar field only in the construction 
of the physical sector of the quantum theory. If we had used that procedure here, the 
factor ordering of the Hamiltonian constraint would have been slightly different. In 
this paper, we started out with the lapse N — \piP2Pi\^^'^ already in the classical theory 
because the resulting quantum Hamiltonian constraint is simpler. In the isotropic case, 
for example, this procedure leads to an analytically soluble model (the one obtained 
in [11] by first starting out with N^i^ = 1, then going to quantum theory, and finally 
making some well-motivated but simplifying assumptions) . It also has some conceptual 
advantages because it avoids the use of "inverse scale factors" altogether. 



D. Simplification of 



grav 



It is straightforward to expand out the Hamiltonian constraint Cgrav using the explicit 
action of operators sin(/ijCj) given by (3.19) and express it as a linear combination of 24 
terms of the type 

Ctj^ VW\Etsgn{Xi) \v\ sgn{Xj)Ef^\, (3.26) 
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(where, i ^ j and as before there is no summation over Unfortunately, the sgn(Ai) 

factors in this expression and in the action of Ef make the result quite complicated. More 
importantly, it is rather difficult to interpret the resulting operator. The expression can be 
simplified if we introduce the volume of V as one of the arguments of the wave function. In 
particular, this would make quantum dynamics easier to compare with that of the Friedmann 
models. With this motivation, let us further re-arrange the configuration variables and set 

v = 2AiA2A3. (3.27) 

The factor of 2 in (3.27) ensures that this v reduces to the v used in the isotropic analysis 
of [10] (if one uses the value of A used there). As the notation suggests, v is directly related 
to the volume of the elementary cell V: 

\> *(Ai, A2, v) = 27r I7I VA |^;| 4i *(Ai, A2, v) . (3.28) 

One's first impulse would be to introduce two other variables in a symmetric fashion, e.g., 
following Misner [53]. Unfortunately, detailed examination shows that they make the con- 
straint (3.20) even less transparent!^ 

Let us simply use Ai,A2,f as the configuration variables in place of Ai, A2, A3. This 
change of variables would be non-trivial in the Schrodinger representation but is completely 
tame here because the norms on H.^^ arc defined using a discrete measure on M^. As a 
consequence, the scalar product is again given by the sum in (3.25), the only difference is 
that A3 is now replaced by v. Since the choice (Ai, A2,f) breaks the permutation symmetry, 
one might have first thought that it would not be appropriate. Somewhat surprisingly, as 
we will now show, it suffices to make the structure of the constraint transparent. (Of course, 
the simplification of the constraint would have persisted if we had chosen to replace either 
Ai or A2 — rather than A3 — with v.) Finally, note that the positive octant is now given by 
Ai > 0, A2 > and v > 0. 

To obtain the explicit action of the constraint, it is extremely convenient to use the fact 
that states ^ in Ti-^^ satisfy the symmetry condition (3.24) and that Cgrav has a well defined 
action on this space. Therefore, to specify its action on any given ^ it suffices to find the 
restriction of the image $(Ai,A2,t') := (Cgrav ^)('^i! ^2; ^^) to the positive octant. The value 
of $ in other octants is determined by its symmetry property. This fact greatly simplifies our 
task because we can use it to eliminate the sgn(Aj) factors in various terms which complicate 
the expression tremendously. 

For concreteness let us focus on one term in the constraint operator (which turns out to 



Misner-like variables — volume and logarithms of metric components — were used in the brief discussion 
of Bianchi I models in [38]. This discussion already recognized that the use of volume as one of the 
arguments of the wave function would lead to simplifications. Dynamics was obtained by starting with 
the Hamiltonian constraint in the /io scheme from [35] and then substituting jji of (3.9) for /^J, in the final 
result. This procedure docs simplify the leading order quantum corrections to dynamics. By contrast, our 
goal is to simplify the full constraint. More importantly, constraint (3.20) is an improvement over that 
of [38] because we introduced pi from the beginning of the quantization procedure and systematically 
defined the operators sm{p.iCi) (in section IIIC). 
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be the most non-trivial one for our simplification): 

(^C2T*)(Ai,A2,v) := (^^E2-sgn(A2) |v| sgn(Ai) ^j"^^*) (Ai,A2,v) 
= [a/H sgn(A2(l )) |'i;-2sgnA2| sgn(Ai) ■\/\v - 2sgnAi - 2sgnA2| ] x 



/t; - 2sgnAi - 2sgnA2 f; - 2sgnA2 \ , . 

* 7, Ai, A2, V - 2sgnAi - sgnA2 . (3.29) 

\ V — 2sgnA2 V / 

If we now restrict the argument of {C^2 ^) the positive octant, the expression simplifies: 
fe-^) = [V^{v-2)^\v-A\] ^(^Ai, ^A2, ^;-4). (3.30) 

\ / + octant V — Z V 

Now the action of this operator is more transparent: the wave function is multiplied by 
functions only of volume and, in the argument of the wave function, volume simply shifts by 
-4 and Ai,A2 are rescaled by multiplicative factors which also depend only on the volume. 
Since the full constraint is a linear combination of terms of this form, its action is also driven 
primarily by volume. As we will see, this key property makes the constraint manageable and 
greatly simplified the task of analyzing the relation between the LQC quantum dynamics 
of Bianchi I and Friedmann Models. From now on, unless otherwise stated, we will restrict 

the argument of the images {Cf^ to lie in the positive octant; its value in other octants is 

given simply by *)(Ai, A2,^;) = ^)(|Ai|, IA2I, |^|). 

The form (3.30) of the action of operators Cf^^ enables us to discuss singularity resolution. 
For completeness, let us first write out the four terms corresponding to i,j=l,2 (which are 
the most complicated of the 24 terms in Cgrav): 

(4++^) (Ai, A2, v) = {v + 2) y^K^^ny ■ ^ C-^K ^A2, + 4) , (3.31) 

(4+-*) (Ai,A2,t;) =t;(^; + 2) •*(^Ai,^A2,^), (3.32) 

(4-+*) (Ai, A2, v) = v{v - 2) • *(^Ai, ^A2, v) , (3.33) 

fe-*) (Ai, A2, v)^(v- 2)V^;|^;-4|- ^(^Ai, ^A2, ^; - 4) . (3.34) 

Recall that, since v is proportional to the volume of the elementary cell, it vanishes when any 
one of the three directional scale factors vanish. Thus, the classical singularity corresponds 
precisely to the points at which v vanishes. Now suppose that the function ^'(Ai, A2, v) has 
no support on points v — ai an initial internal time Tg. As it evolves via (3.22), can it end 

up having support on such points? We will argue that this is impossible. 

Let us decompose 7i§^^ as 7i§„^ = TYfj^^J © '^rcg^ where \l/(Ai,A2,f) is in Tifj^g if it 
has support only on points with v = and it is in Tifl^^ if it has no support on points 

with V — 0. Now, all the operators C^^ have a factor of ^/v acting on the right (see Eq. 

(3.26)). It ensures that each C^^ annihilates every state in Hfl^^. Therefore Hfl^^ is left 
invariant by the evolution. More importantly, because of the pre-factors oi v ±2 and v ±4 
the action of the 4 operators in (3.31) - (3.34) preserves H%g^. This property is shared 
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also by Cf^^ for other values of i,j and hence by Cgrav and all its powers.^ Therefore, the 



relational dynamics of (3.22) decouples TY^j^^g from Tif^g^. In particular, if one starts out 
with a "regular" quantum state at T = 0, it remains regular throughout the evolution. In 
this precise sense, the singularity is resolved. 

Next, let us write out explicitly the full Hamiltonian constraint (3.22): 

*(Ai, A2, v; T) =^v^[(^; + 2) *+(Ai, A2, v- T) - (v + 2) *+(Ai, A2, v; T) 

-{v- 2) ^0 (Ai, A2, V- T) + {v- 2)Vl^^-4|^4(Ai, A2, v- T) 

'(3.35) 

where ^^^4 are defined as follows: 

-1-4 7; + 2 \ fv + A 

*±(Ai,A2,^;r)=* — --Ai, A2,^±4;r — -•Ai,A2,^±4;r 

\v ± 2 V J \v ±2 

+ *(— •A.^^-A„.±4;r)+*(!i±l.A.,A„.±4;r 

\ V V ±2 J \ V 

f V ±2 \ f V ±A \ 

+ ^ Ai, A2, w ± 4; T + ^ Ai, —— ■ A2, V ± 4; T , (3.36) 



v±2 



and 



*±(Ai, A2, T) = * (^-^ ■ Ai, ^ • A2, v-T^+m {^-^ ■ Ai, A2, v; T 

+ * {-^ • Ai, ^ • A2, ^; r) + * (-^ ■ Ai, A2, V- T 
\v±2 V J \v ±2 

+ ^ (ai, • A2, v-T^+^> (^Ai, ^ • A2, V- , (3.37) 

where, as before, we have given the restriction of the image of Cgrav to the positive octant. 
Because Hf^^ is left invariant by evolution we can in fact restrict Ai,A2,v to be strictly 
positive. On the right sides of (3.36) and (3.37), arguments of ^ can take negative values. 
However, since ^(Ai, A2,f) = ^'(lAil, IA2I, |f|), we can just introduce absolute value signs on 
these arguments. Consequently, knowing the restriction of \t' to the positive octant, (3.36) 
and (3.37) enable us to directly calculate its image under Cgrav In particular, numerical 
evolutions can be carried out by restricting oneself to the positive octant. 

Let us now examine the structure of this equation. As in the isotropic case, the right side 
is a difference equation. As far as the v dependence is concerned, the steps are uniform: the 
argument of the wave function involves v — 4, v + 4 exactly as in the isotropic case. The 
step sizes are also the same as in [10] because, as noted above, our variable v is in precise 
agreement with that used in the isotropic case. There is again superselection. For each 
e e [0, 4), let us introduce a 'lattice' >Cg consisting of points v = 4n if e = and v — 2n-\- e 



^ To make this argument mathematically rigorous one would have to establish that Cgrav is essentially 
self-adjoint and its self adjoint extension also shares this property (or a suitable generalization thereof). 
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if e 7^ 0.^ Then the quantum evolution — as well as the action of the Dirac obscrvablcs — 
preserves the subspaces "H^^^ consisting of states with w-support on /Z^. The most interesting 
of these sectors is the one labelled by e = since it contains the classically singular points, 
V — 0. Therefore in what follows, unless otherwise stated, we will restrict ourselves to this 
sector. 

The dependence of Cgrav ^ on Ai, A2, by contrast, is much more difficult to control tech- 
nically because the first two arguments of the wave function cannot be chosen to lie on a 
regular lattice in any simple way. In particular, even if we started out with a wave function 
which has support only on a lattice, say Ai = nXg for some Aq, the action of Cgrav shifts 
support to points such as Ai = [{v±2)/v]nXo which do not lie on this lattice. Thus, there is 
no obvious superselection with respect to Ai and A2; we have to work with the entire they 
span. Had it been permissible to set /Xj oc / 1/ |pi | , we could have restricted A^ to lie on a reg- 
ular lattice [36] . Then, following [40] , we could have repeated the strategy used successfully 
in the isotropic case in [11] to simplify dynamics by carrying out a Fourier transform to pass 
to variables which are conjugate to Ai, A2. However, as remarked earlier, that choice of /Xj is 
inadmissible and hence the strategy cannot be repeated in the Bianchi I case. Nonetheless, 
it is still feasible to carry out numerical simulations. For, if one knows the support of the 
quantum state at an initial time Tq and the number of time-steps across which one wants 
to evolve, one can calculate the number of points on a (irregular) grid in the A1-A2 plane on 
which the wave function will have support. Numerical work has in fact already commenced 
[41]. It would be interesting to investigate whether the efficient algorithms that have been 
introduced in the context of regular lattices [54] can be extended to this case. 

We will conclude this discussion by noting that it is possible to read off some qualitative 
features of dynamics from (3.35) - (3.37). Since the steps in v of this difference equation 
are the same as those in the isotropic case, the dynamics of volume — and also of the matter 
density Pmatt, since p(T) is a constant of motion — would be qualitatively similar to that in 
the isotropic case. What about anisotropies? The Xj {I — 1,2) do not feature in the overall 
numerical factors in (3.35); they appear only in the argument of the wave functions. Under 
the action of Cgrav, these arguments get rescaled by factors v ±4/v ±2, v±2/v and v/v±2. 
For large volumes, or more precisely low densities, these factors go as 1 -|- 0(pmatt/ppi)- 
Hence, to leading order, we will recover of the classical result that a\a2a^{Hi — Hj) are 
constants, where are the directional scale factors and Hi :— dlncj/dt, the directional 
Hubble parameters. Since quantum corrections go as p/ ppi they are utterly negligible away 
from the Planck regime. 

In the next section we will discuss three important features of dynamics dictated by (3.35) 
which provide significant physical intuition in complementary directions. 

IV. PROPERTIES OF THE LQC QUANTUM DYNAMICS 

This section is divided into three parts. Since we have used the same general procedure 
as in the isotropic case it is natural to ask how the quantum dynamics of (3.35) compares 
to that in [10]. In the first part we show that there is a natural projection from a dense 
subspace of the physical Hilbert space of the Bianchi I model to that of the Friedmann model 



^ As in the isotropic case, the lattice is doubled if e ^ or 2 because of the symmetry property of our wave 
functions. 
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which maps the Bianchi I Hamiltonian constraint to that of the Friedmann modeL This 
result boosts confidence in the overall coherence and reliability of the quantization scheme 
used in LQC. In various isotropic models [10, 14-16, 18], one can derive certain effective 
equations. Somewhat surprisingly, for states which are semi-classical at a late initial time, 
they faithfully capture quantum dynamics throughout the entire evolution, including the 
bounce. The same considerations lead to effective equations in Bianchi I models which were 
already analyzed by Chiou and Vandersloot in Appendix C of [37] . In the second sub-section 
we briefly discuss these equations and their consequences. In the third, we show that, as 
in the isotropic case [10, 11], there is a precise sense in which the LQC quantum dynamics 
reduces to that of the Wheeler- De Witt theory in the low curvature regime. 



A. Relation to the LQC Friedmann dynamics 

The problem of comparing dynamics of a more general system with that of a restricted, 
symmetry reduced one has been discussed in the literature in several contexts. In the classical 
theory, symmetric states often provide symplectic sub-manifolds F^es of the more general 
phase spaces FGcn- Furthermore Frcs are preserved by the dynamics on Fccn- Therefore, it is 
tempting to repeat the same strategy in the qiiantum theory. Indeed, sometimes it is possible 
to find natural sub-spaces of states with additional symmetry in the full Hilbert space 
T^Gen of the more general system. However, generically T^Res is not left invariant by the more 
general dynamics (see, e.g., [43, 44]). In our case, one can introduce an isotropic sub-space 
of T^Res in the quantum theory based on any given fiducial cell V: isotropic states correspond 
to wave functions \I'(Ai, A2, which have support only at points Ai = A2 = (f/2)^/^. (But 
note that this sub-space is not invariantly defined; it is tied to V!) It is easy to check that 
the space Hrcs of these states is not left invariant by the Bianchi I quantum dynamics (3.35). 

However, this fact cannot be interpreted as saying that there is no simple relation between 
the quantum dynamics of the two theories: since restriction to Tt^^s amounts to a sharp 
freezing of anisotropic degrees of freedom, in view of the quantum uncertainty principle, 
this procedure is not well suited to compare the quantum dynamics of the two systems. As 
pointed out in section I, a better strategy is to integrate out the extra, anisotropic degrees 
of freedom. This would correspond to a projection map from Ticen to T^Res rather than an 
embedding of T^Res into Ticen- 

Consider first, as an elementary example, a particle moving in R^. Suppose that the 
potential depends only on z so that dynamics has a symmetry in the x, y directions. In 
the classical theory, there are several natural embeddings of the phase space T^es 

into FQen- 

For example, we can set {z,Pz) —>■ {x=Xo, y=yo,z; Px=0, Py=0, Pz) and the Hamiltonian 
vector field of the full theory is then tangential to the images of each of these embeddings. 
However, in the quantum theory the Hilbert space T^oen of the full system is L^(R^,d^x) 
and there is no natural embedding %l){z) '^{x^y^z). The classical strategy would suggest 
setting \&(a;, y, z) = 6{x, Xo) S{y, yo)ip{z) but this is not a normalizable state in Ticcn for any 
ip{z). Even if one were to ignore this fact and try to evolve these states, one would find that 
they are not preserved by the full Hamiltonian operator H. 

Note however that there is a natural projection P from a dense subspace in TYcen to that 
in T^Res: 

^{x,y,z) ^ {m){z) -.^ [dx [ dy^{x,y,z) = jP{z) . (4.1) 
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(For example, we can choose the dense subspace to be the space of smooth functions of 
compact support.) Furthermore, under this projection, the Hamiltonian operator 

H = -{h^/2m)A + V{z) 

of the general system is mapped to the Hamiltonian operator 

h:^-{hy2m)dydz'' + V{z) 

of the reduced system. Hence solutions \E'(a;, t) of the Schrodinger equation of the full system 
are mapped to solutions ip{z,t) of the reduced system. Finally, this projection strategy 
continues to work for more general Hamiltonians of the type P{z)pi + V{z) which again 
have a symmetry in the x, y directions. 

Let us return to the Bianchi I model and define a projection P from states ^'(Ai, A2, v) of 
the Bianchi I model to the states ijj{v) of the Friedmann model of [10] as follows: 

*(Ai, A2, v) ^ (P*)(^;) := ^2, v) = V'(^) ■ (4-2) 

Ai,A2 

(The idea of using such a map already appeared in [42] where the map was defined between 
elements of Cyl* of the locally rotationally symmetric Bianchi 1 model and that of the 
Friedmann model.) Again, P is a well defined projection from a dense subspace of the 
Bianchi I Hilbert space to a dense subspace of the Friedmann Hilbert space, consisting, for 
example, of states which have support only on a finite number of points. As is manifest 
from (4.2), its effect is to focus on volume by "integrating out" the anisotropic degrees of 
freedom with the same volume. Applying this projection map P to Eq. (3.35), we find 

dlip{v-T) = 37iG \{v + 2)^/v{v + 4)^(u + 4; T) - 2v^ij{v; T) 



+ (v-2)y^;<^7^^V'(v-4;r) . (4.3) 

This is precisely the quantum constraint describing the LQC dynamics of the Friedmann 
model with lapse^ N — \pf'^'^ ■ The reason for the exact agreement is two-fold. First, the 
Hamiltonian constraint Cgrav of the Bianchi I model is a difference operator whose coefficients 
depend only on v and, second, the shift in the argument is dictated only by v. Thus, 
conceptually, Ai, A2 are "inert directions" in the same sense that x,y are in the elementary 
example discussed above. To summarize, there is a simple — and exact — relation between 
quantum dynamics of the two theories. It would be interesting to investigate if this result 
admits a suitable extension to other Bianchi models [33, 55]. 

In completely general situations, of course, this exact agreement will not persist: the 
projected dynamics will provide extremely non-trivial corrections to the dynamics of the 
simpler system. However, the BKL conjecture says that the dynamics of general relativity 
greatly simplifies near space-like singularities: In this regime, the time evolution at any one 
spatial point is well modelled by that of Bianchi I cosmology. Therefore, in a large class of 



^ As noted at the end of section III C, the analysis in [10] began with the lapse = 1 and therefore leads to 
a slightly different factor ordering. Had one used N = from the beginning as in the current paper, 

one would have obtained the factor ordering used in [11]. Eq. (4.3) matches exactly with that constraint. 
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situations there may well be a sense in which the quantum dynamics in the deep Planck 
regime can be projected to that of the Friedmann model with only small corrections. If so, 
the Planck scale quantum dynamics of the isotropic, homogeneous degree of freedom in the 
full theory will be much simpler than what one would have a priori expected. 

B. Effective equations 

Physically, the most interesting quantum states are those that are sharply peaked at 
a classical trajectory at late times. As explained in section I, in the isotropic case such 
states remain peaked at certain effective trajectories at all times, including the epoch during 
which the universe undergoes a quantum bounce. Thus, even in the deep Planck regime 
quantum physics is well captured by a smooth metric although its dynamics can no longer 
be approximated by the classical Einstein's equations and its components now contain large, 
/i-dependent terms. The effective equations obeyed by these geometries were first derived 
using ideas from geometrical quantum mechanics [56, 57]. However, the assumptions made 
in these derivations break down in the deep Planck regime. Therefore a priori there was no 
reason to expect these equations to describe quantum dynamics so well also in the Planck 
regime. That they do was first shown by numerical simulations of the exact quantum 
equations [9, 10] in the k=0, A=0 case. It was then realized that this model is in fact exactly 
soluble [11, 51] and the power of the effective equations could be attributed to this property. 
However, k=0 models with non-zero cosmological constant and the closed k=l models do 
not appear to be exactly soluble. Yet, numerical solutions of the exact quantum equations 
show that the effective equations continue to capture full quantum dynamics extremely well 
[14-16]. 

New light was shed on this phenomenon by recent work on a path integral formulation 
of quantum cosmology [58]. The idea here is to return to the original derivation of path 
integrals due to Feynman and Hibbs [59] starting from quantum mechanics. In the isotropic 
case, then, the strategy is to begin with the kinematics and dynamics of LQC and then 
rewrite the transition amplitudes as path integrals. The resulting framework has several 
novel features. First, because the LQC kinematics relies on quantum geometry, paths that 
feature in the final integral are different from what one would have naively expected from 
the Wheeler-DeWitt theory. Second, the action that features in the measure is not the 
Einstein-Hilbert action but contains non-trivial quantum corrections. When expressed in 
the phase space language, L — pq — H[p,q), the "Hamiltonian" H turns out to be precisely 
the effective Hamiltonian constraint derived in [56, 57], even though this casting of the LQC 
transition amplitudes in the path integral language is exact and does not pre-suppose that we 
are away from the Planck regime. Now, in the path integral approach, we have the following 
general paradigm. Consider the equations obtained by varying the action that appears in 
the path integral. (Generally these are just the classical equations but in LQC they turn 
out to be the effective equations of [10, 16, 57].) Fix a path representing a solution to these 
equations. If the action evaluated along this path is large compared to h then that solution 
is a good approximation to full quantum dynamics. If one applies this idea to isotropic 
LQC, one is led to conclude that solutions to the effective equations of [56, 57] should be 
good approximations to full quantum dynamics also in the k=0, A ^0 and k=l cases. This 
is precisely what one finds in numerical simulations. Thus, the path integral approach may 
well provide a deeper explanation of the power of effective equations. While such a path 
integral analysis is yet to be carried out in detail in the anisotropic case, because of the 
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situation in the simpler cases it is of considerable interest to find effective equations and 
study their implications. 

This task was carried out already by Chiou and Vandersloot in the Appendix C of [37]. 
We will summarize the relevant results and briefiy comment on the general picture that 
emerges. 

Without loss of generality, we can restrict ourselves to the positive octant. Then the 
effective Hamiltonian constraint is given simply by the direct classical analog of (3.14): 

P(T)+Cv = (4-4) 

where 

C?av = ~ g7rG'7^A ^ ^™ ^^^^ '"^^^ ^ ^™ ^'^^'^ ^™ ^^^^ ^ ^^^^ ^™ ^^^^^ ' 

Since sin a; is bounded by 1 for all ,x, these equations immediately imply that the matter 
density, Pmatt = P(t)/2V"^ = P(T)/2piP2P3 can never become greater than the critical density 
Pcrit ~ 0.41pp], first found in the isotropic case [10-12, 16, 18]. Since p becomes infinite 
at the big bang singularity in the classical evolution, there is a precise sense in which the 
singularity is resolved in the effective theory. 

Effective equations are obtained via Poisson brackets as in section II but using (4.4) in 
place of the classical Hamiltonian constraint. This gives, for example, 

dpi PWP1P2P3 , . _ \ . . 

— = — — cos(/xiCi) sm//2C2 + sm/i3C3 , (4.6) 



and 

dci P2P3 



dr A-fi- 



sm piCi sm H2C2 + sm piCi sm p^c^ + sm /i2C2 sm p^c^ 
H -— cos niCi [ sm /X2C2 + sm H3C3 ) — cos /i2C2 ( sm niCi + sm /isCs j 



cos//3C3(sin//iCi + sin/X2C2) . (4.7) 

Equations for p2, C2 and ps, C3 are obtained by cyclic permutations. These effective equations 
include "leading order quantum corrections" to the classical evolution equations (2.13) and 
(2.14). In any solution, these corrections become negligible in the distant past and in the 
distant future. As we noted in section II, the shear S defined in Eq. (2.21) is a constant of 
motion in the classical theory. This is no longer the case in the effective theory. However, one 
can show that it remains finite throughout the evolution, becomes approximately constant 
in the low curvature region both in the distant past and in the distant future. Furthermore, 
its value in the distant future is the same as that in the distant past along any effective 
trajectory in the phase space. 

Vandersloot (personal communication) has also carried out numerical integration of these 
equations. In the isotropic case each effective trajectory undergoes a quantum bounce when 
the matter density Pmatt achieves a critical value pcrit ~ 0.41ppi. As one might expect, 
now the situation is more complicated because of the additional degrees of freedom. First, 
there are now several distinct "bounces" . More precisely, in addition to Pmatt {or the scalar 
curvature), we now have to keep track of the three Hubble rates Hi which directly control 
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the Weyl curvature. In the backward evolution towards the classical big bang, Einstein's 
equations approximate the effective equations extremely well until the density of one of the 
Hi enters the Planck regime. Then the quantum corrections start rising quickly. Their 
net effect is to dilute the quantity in question. Once the quantity exits the Planck regime 
as a result of this dilution, quantum geometry effects again become negligible. Thus, as 
in the isotropic case, one avoids the ultraviolet-infrared tension [21] because the quantum 
geometry effects are extremely strong in the Planck regime but die off extremely quickly as 
the system exits this regime. Secondly, the "volume" or the "density bounce" occurs when 
the matter density is lower than pcrit- This is not surprising because what matters is the 
total energy density and now there is also a contribution from gravitational waves. Finally, 
although there are distinct "bounces" for density (or scalar curvature) and the Hi (or the 
Weyl curvature invariants), they all occur near each other in the relational time T. 

There are indications that the general scenario provided by effective equations correctly 
captures the qualitative features of the full quantum evolution. However, the arguments are 
not conclusive. For conclusive evidence for (or against) this picture, one needs numerical 
simulations [41] of the exact quantum equations of section III D, or a detailed, path integral 
treatment of the Bianchi I models along the lines of [58] . 



C. Relation to the Wheeler-DeWitt Dynamics 

Quantum dynamics of LQC is governed by a difference — rather than a differential — 
equation because of the quantum geometry effects. However, we will now show that, as 
in the isotropic case [10, 11, 16], the LQC quantum dynamics is well approximated by the 
Wheeler-DeWitt (WDW) differential equation away from the Planck regime where quantum 
geometry effects become negligible. 

In the WDW theory the directional scale factors, and hence the three Aj can assume any 
real value and it is simpler to work with the three A.; rather than with Ai, A2, v = 2A1A2A3. 
Let us therefore set ^^.(Ai, A2, A3; T) = ^^(Ai, \2:V; T) and assume that ^ admits a smooth 
extension to all real values of Aj. The idea is to pair various terms in Eqs. (3.36) and (3.37) 
in such a way so that two of the three arguments of ^ are the same. For example, one such 
pair is 

^f^-Ai,^-A2,A3;T') and ^ • A^, ^ • A^, A3; tV (4.8) 



v + 2 V - y — \^^ + 2 

Next, let us define v' — v + 2 and A2 = v'X2/{v' — 2) so that we have 




aA7+4 = y^7T2 = J Ai + -— A'2A3 . (4.9) 



Ignoring the common pre-factors in Eqs. (3.36) - (3.37), the two paired terms in Eq. (4.8) 
can be expressed as: 

^AyT2 ^(Ai + j^, A'2, A3; T) - V^T^ ^(Ai - 3;^, A'2, A3; T) 
= afr *(Ai, A'2, A3; T) + o((^)- V^*) 
= ^ ^v^3^(Ai, K, A3; T) + 0((^)" 4 v^*) (4.10) 
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where n > 1. (Notice that the v' in the denominator in front of the partial derivative will 
cancel the v + 2 pre-factor in Eq. (3.35).) One can suitably pair all terms in (3.36) and 
(3.37) and express them as differential operators with corrections which are small for large 
values of Aj. Let us ignore these corrections — i.e. assume that the {1/ XiXj)"'d^ ^/v'^ is 
negligible for n > 1 because ^ is slowly varying and we are in the low density, large scale- 
factor regime. Then we find that the LQC Hamiltonian constraint (3.35) reduces to a rather 
simple differential equation: 



9^mAi,A2,A3;T)=^ 



d d d d d d d d 

Ai7TT~A27rT h Ai A3— h A2 ^, Ai— h A2 ^, A3 



dXi 8X2 dXi 8X3 8X2 8X1 8X2 8X3 
This equation can be further simplified by introducing ai — log A^ and $ = V^^- The result 



is: 



f)2 f)2 o2 

4$(ai, a2, as; T) = 167rG [— -- + — -- + -—-1 $(ai, ^2, a,; T) , (4.12) 

i8ai8a2 8ai8a3 8a28a2,^ 

where v is now given by 2 exp(^ ai). This is precisely the equation we would have obtained if 
we had started from the classical Hamiltonian constraint, used the Schrodinger quantization 
and the "covariant factor ordering" of the constraint as in the WDW theory. Thus, the 
LQC Hamiltonian constraint reduces to the WDW equation under the assumption that 
is slowly varying in the sense that {l/XiXjY8^^^ can be neglected for n > 1 relative to 
the term for n = 1. Since (AjA^)^ is essentially the area of the i-j face of the fiducial cell 
V in Planck units, this should be an excellent approximation well away from the Planck 
regime. However, in the Planck regime itself the terms which are neglected in the LQC 
dynamics are comparable to the terms which are kept whence, as in the isotropic case, the 
WDW evolution completely fails to approximate the LQC dynamics. 



V. DISCUSSION 

In this paper we extended the "improved" LQC dynamics of Friedmann space-times [10] 
to obtain a coherent quantum theory of Bianchi I models. As in the isotropic case, we 
restricted the matter source to be a massless scalar field since it serves as a viable relational 
time parameter (a la Leibniz) both in the classical and quantum theories. However, it is 
rather straightforward to accommodate additional matter fields in this framework. 

To incorporate the Bianchi I model, we had to overcome several significant obstacles. 
First, using discrete symmetries we showed that to specify dynamics it suffices to focus 
just on the positive octant. This simplified our task considerably. Second, in section HI B 
we introduced a more precise correspondence between LQG and LQC and used it to fix the 
parameters jli that determine the elementary plaquettes, holonomies around which define the 
curvature operator Fab'' ■ This procedure led us to the expressions fil = ( |pi | A £pj) / \p2P3 \ , etc. 
They reduce to the expression jji^ — (A£|j)/|p| of the isotropic models [10, 16, 18]. But even 
there, the current reasoning has the advantage that it uses only quantum geometry, avoiding 
reference to classical areas even in the intermediate steps. However, because of this rather 
complicated dependence of /ij on pi, the task of defining operators sin/ijQ seems hopelessly 
difficult at first. Indeed, this was the key reason why the earher treatments [36, 37, 40] took 
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a short cut and simply set p,'^ = {Aip^)/\pi\ by appealing to the relation /i^ = {A£p^)/\p\ 
in the isotropic case. With this choice, quantization of the Hamiltonian constraint became 
straightforward and the final Bianchi I quantum theory resembled three copies of that of 
the Priedmann model. However, this result had the physically unacceptable consequence 
that significant departures from general relativity could occur in "tame" situations. By a 
non-trivial extension of the geometrical reasoning used in the isotropic case, in section 111 C 
we were able to define the operators sin /ijCj for our expressions of jli. However, the structure 
of the resulting Hamiltonian constraint turned out to be rather opaque. To simplify its form, 
in section HI D we introduced volume as one of the arguments of the wave functions. The 
action of the gravitational part of the Hamiltonian constraint then became transparent: it 
turned out to be a difference operator where the multiplicative coefficients in individual 
terms depend only on volume and the change in the arguments of the wave functions also 
depends only on volume; individual anisotropics do not feature (see (3.35) - (3.37)). This 
simplification enabled us to show that the sector Hfl^ of quantum states which have no 
support on classically singular configurations is preserved by quantum dynamics. In this 
precise sense the big-bang singularity is resolved. Furthermore, this quantum dynamics is 
free from the physical drawbacks of the older scheme mentioned above. 

In section IV we explored three consequences of quantum dynamics in some detail. First, 
wc showed that there is a projection map P : Ticcn '^Rcs from the Hilbcrt space of the 
more general Bianchi 1 model to that of the more restricted Friedmann model which maps 
the Bianchi I quantum constraint exactly to the Friedmann quantum constraint. This is 
possible because, as noted above, it is just the volume — rather than the anisotropies — 
that govern the action of the Bianchi I quantum constraint. This result is of considerable 
interest because, in view of the BKL conjecture, it suggests that near generic space-like 
singularities the LQC of Friedmann models may capture qualitative features of the full, 
LQG dynamics of the isotropic, homogeneous degree of freedom. In section IV B we briefly 
recalled the effective equations of Chiou and Vandersloot (see Appendix C of [37]). These 
equations provide intuition for the rich structure of quantum bounces in the Bianchi I model. 
Their analysis suggests that classical general relativity is an excellent approximation away 
from the Planck regime. However, in the Planck regime quantum geometry effects rise 
steeply and forcefully counter the tendency of the classical equations to drive the matter 
density, the Ricci scalar and Weyl invariants to infinity. (In particular, as in the isotropic 
case, the matter density is again bounded above by pcrit ~ 0.41ppi.) Thus the quantum 
geometry effects dilute these quantities and, once the quantity exits the Planck regime, 
classical general relativity again becomes an excellent approximation. In section IV C we 
showed that, as in the isotropic case [10, 11, 16], there is a precise sense in which LQC 
dynamics is well approximated by that of the WDW theory once quantum geometry effects 
become negligible. 

The rather complicated dependence of fii on is also necessary to remove a fundamental 
conceptual hmitation of the older treatments of the Bianchi I model. Recall that, because we 
have homogeneity and the spatial topology is non-compact, we have to introduce a fiducial 
cell V to construct a Lagrangian or a Hamiltonian framework. Of course, the final physical 
results must be independent of this choice. At first this seems like an innocuous requirement 
but it turns out to be rather powerful. We will now recall from [40] the argument that this 
condition is violated with the simpler choice /if = {A£p^)/\pi\ but respected by the more 
complicated choice we were led to from LQG. 

For definiteness, let us fix a fiducial metric °qab and denote by Lj the lengths of the 
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edges of the fiducial cell V. Suppose wc were to use a different cell, V whose edges have 
lengths L\ = j3iLi (no summation over i). Since the basic canonical fields and are 
insensitive to the choice of the cell, Eq. (2.3) implies that the labels q and pi we used to 
characterize them change to c'^ = /3iCi, p'l = P2P3P1, etc. The gravitational part of the 
classical Hamiltonian constraint (2.12) is just rescaled by an overall factor and 
the inverse symplectic structure is rescaled by (/9i/52/33)~^. Hence the Hamiltonian vector 
field is rescaled by exactly as it should because the lapse is rescaled by the same 

factor. Thus, as one would expect, the classical Hamiltonian flow is insensitive to the 
change V — > V. What is the situation in the quantum theory? Physical states belong to 
the kernel of the Hamiltonian constraint operator Ch whence the two quantum theories will 
carry the same physics only if Ch is changed at most by an overall rescaling. Analysis is a 
bit more involved than in the classical case because Cgrav involves factors of sin/ijQ. Now, 
under V V, our p,i transform as p,i p'l = Pi^fli, whence fi'ic[ = fiiCi, etc, and the 
Hamiltonian constraint (3.14) is rescaled by an overall multiplicative factor {l3il32l3'iY just 
as in the classical theory. What happens if we set pL^ = A£pj/|pj| as in [36, 37, 40]? Then, 
we arc led to fL[c[ = (/5i/V/^2/93) fiiCi etc. Since the constraint (3.14) is a sum of terms 
of the type P1P2IP3I sin/iiCi sin/i2C2 it has a rather uncontrolled transformation property 
and is not simply rescaled by an overall factor. It is then not surprising that, in the Planck 
regime, the dynamical predictions of the resulting quantum theory (as well as of the effective 
theory) depend on the choice of the elementary cell. It is rather remarkable that the more 
complicated form of jli that we are led to from LQG kinematics has exactly the right form 
to make quantum dynamics insensitive to the choice of the fiducial cell V. As mentioned 
above, it also ensures that the predictions of quantum theory is free of drawbacks of the 
earlier treatments [36] , such as the correlation between the bounce and "directional densities" 
which do not have an invariant significance. 

Prom physical considerations, as in the isotropic case, it would be most interesting to 
start at a "late time" with states that are sharply peaked at a classical solution in which the 
three scale factors assume values for which the curvature is "tame" and pi^r) is very large 
compared to h in classical units c=G=l. One would then evolve these states backward and 
forward in the "internal" time T. As we just discussed, analytical considerations show that, 
since the initial wave function is in Hfl^, it will continue to be in that sub-space; there is 
no danger that the expectation values of curvature, anisotropics or density would diverge. 
But several important questions remain. Are there quantum bounces with a pre-big-bang 
branch again corresponding to a large, classical universe in the distant past? Is there is a 
clear distinction between evolutions of data in which there are significant initial anisotropics 
and data which represent only perturbations on isotropic situations? Even in the second 
case, do anisotropics grow (or decay) following predictions of the classical theory or are 
there noticeable deviations because of accumulations of quantum effects over large time 
periods? Numerical simulations of the LQC equations are essential to provide confidence in 
the general scenario suggested by effective equations and to supply us with detailed Planck 
scale physics. 

Finally, let us return to full LQG. At the present stage of development, there appears 
to be considerable freedom in the definition of the quantum Hamiltonian constraint in the 
full theory. Furthermore, our current understanding of the physical implications of these 
choices is quite limited. Already in the isotropic models, the "improved" dynamics scheme 
provided some useful lessons: it brought out the fact that these choices can be non-trivially 
narrowed down by carefully analyzing conceptual issues (e.g., requiring that the physical 
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results should be independent of auxiliary structures introduced in the intermediate steps) 
and by working out the physical consequences of the theory in detail (to ensure that the 
quantum geometry effects are not dominant in the low energy regime). Rather innocuous 
choices — such as those made in arriving at the older "/^o-scheme" — can lead to unaccept- 
able consequences on both these fronts [12]. The Bianchi I analysis has sharpened these 
lessons considerably. The fact that the kinematical interplay between LQG and LQC has a 
deep impact on the viability of quantum dynamics is especially revealing. A quantum anal- 
ysis of inhomogeneous perturbations around Bianchi I backgrounds is therefore a promising 
direction for understanding the physical implications of the choices that have to be made 
in the definition of the Hamiltonian constraint in full LQG. Such an analysis is likely to 
narrow down choices and lead us to viable quantization schemes in LQG that lead to a good 
semi-classical behavior. 
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APPENDIX A: PARITY SYMMETRIES 

In this appendix we recall and extend results on parity symmetries obtained in [46] . 

In non-gravitational physics, parity transformations are normally taken to be discrete 
diffeomorphisms the physical space which are isometrics of the flat 3-metric 

thereon. In the phase space formulation of general relativity, we do not have a flat metric — 
or indeed, any fixed metric. However, if the dynamical variables have internal indices — such 
as the triads and connections used in LQG — we can use the fact that the internal space / 
is a vector space equipped with a fiat metric qij to define parity operations on the internal 
indices. Associated with any unit internal vector C,^ , there is a parity operator H^ which 
reflects the internal vectors in the 2-plane orthogonal to ^. This operation induces a natural 
action on triads e", the connections and the conjugate momenta =: (l/87rG7)£'f (since 
they are internal vectors or co- vectors). It turns out that e" are proper internal co- vectors 
while and J'f arc pseudo internal vectors and co-vectors, respectively. These geometrical 
considerations show that the Barbcro-Immirzi parameter 7 must change sign under any one 
of these parity operations, i.e., if it has the value I7I for say, positively oriented triads, it 
should have the value — I7I for negatively oriented triads. Its value on degenerate triads is 
ambiguous so on the degenerate sector we cannot unambiguously recover the triads from 
the momenta P,"-. If one were to make 7 a dynamical field [61, 62], it follows that the field 
should be a pseudo-scala.i under internal parity transformations; geometrical considerations 
involving torsion have led to the same conclusion in [62]. (For details, see [60]). 

In the diagonal Bianchi I model, we can restrict ourselves just to three parity operations 
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Hj. Under their action, the canonical variables Cj,pj transform as follows: 



ni(Ci, C2, C3) = (Ci, -C2, -C3) 



ni(pi,P2,P3) 



-Pl,P2,P3) 



(Al) 



and the action of 112,113 is given by cyclic permutations. Under any of these maps Ilj, 
the Hamiltonian (2.12) is left invariant. This is just as one would expect because Ilj are 
simply large gauge transformations of the theory under which the physical metric Qab and 
the extrinsic curvature Kab do not change. It is clear from the action (Al) that if one 
knows the dynamical trajectories on the octant Pi > oi the phase space, then dynamical 
trajectories on any other octant can be obtained just by applying a suitable (combination 
of) Ilj. Therefore, in the classical theory one can restrict one's attention just to the positive 
octant. 

Let us now turn to the quantum theory. We now have three operators Ilj. Their action 
on states is given by 

ni^(Ai, A2, A3) = ^(-Ai, A2A3) , (A2) 
etc. What is the induced action on operators? Since 



niAini^(Ai, A2, A3) = Ui (^Ai ^(-Ai, A2, A3) 

= -Ai^'(Ai,A2,A3), 



(A3) 



we have 



HiAiHi = -Al. 



(A4) 



The Hamiltonian constraint operator is given by Eqs. (3.20)- (3.21). To calculate its trans- 
formation property under parity maps, in addition to (A4), we also need the transformation 
property of operators sin /ijQ. An inspection of Eq. (3.21) shows that, in view of the Bianchi 
I symmetries, it is sufficient to calculate itj sin/iicilli. We have: 

U^ sin ^l^c,U,^{X^, A2, A3) ^ ^Ui h{-X, - ^-^^^tM^ A2, A3) - *(-Ai + ^l^^tM^ A2, A3)' 

Zl L A2A3 A2A3 



1 r 

2i 



*(^^ - - *(Ai + A2, A3) 

A2A3 A2A3 



sgn(Ai 



= sin/iiCi^'(Ai, A2, A3), 



whence 



ill sin /xiCiIIi = sin /xiCi. 
An identical calculation shows that 



(A5) 

(A6) 



n2sin/xicin2*(Ai, A2, A3) = ^n2 



*(^^-7^'-^^'^3)-*(Ai + 
A2JA3 



sgn(Ai) 

(-A2)A3 



— A2, A3) 



1 

2i 



*(A, + ^.A..A,)-*(A,-«^.A,.A3) 



A2A3 

-sin/iiCi*(Ai, A2, A3) , 



A,A 



2^3 



(A7) 



and similarly for 113. Therefore, we have: 

II2 sin p,iCiTl2 — — sin fiiCi, and Il3sin fiiCitls — — sin fliCi. 



(A8) 
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These transformation properties of sin /iiCi under Ilj simply mirror the transformation prop- 
erties of Ci under the three parity operations Hi in the classical theory. (Note that, because 
of the absolute value signs in the expressions (3.9), p,i do not change under any of the parity 
maps.) 

Prom Eqs. (3.20)- (3.21) it now immediately follows that the gravitational part of the 
Hamiltonian constraint is left invariant under flj. Since p^.^ is manifestly invariant, we have: 

Ch fli = Ch. (A9) 

just as in the classical theory. Because of this invariance property, given any state ^ e "^lin^; 
the restriction to the positive octant of its image under Cgrav determines its image everywhere 
on Ti^in^- A^ section HID, this property simplifies the task of finding the explicit 

action of the Hamiltonian constraint considerably. 
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